d-wave Superconductivity in the Hubbard Model 
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The superconducting instabilities of the doped repulsive 
2D Hubbard model are studied in the intermediate to strong 
coupling regime with help of the Dynamical Cluster Approxi- 
mation (DC A). To solve the effective cluster problem we em- 
ploy an extended Non Crossing Approximation (NCA) , which 
allows for a transition to the broken symmetry state. At suffi- 
ciently low temperatures we find stable d-wave solutions with 
off-diagonal long range order. The maximal T c ~ 150K oc- 
curs for a doping 8 « 20% and the doping dependence of the 
transition temperatures agrees well with the generic high-T c 
phase diagram. 



Introduction The discovery of high-T c superconduc- 
tors has stimulated strong experimental and theoretical 
interest in the field of strongly correlated electron sys- 
tems. After a decade of intensive studies we are still 
far from a complete understanding of the rich physics 
observed in high-T c cuprates [Ej. Angle resolved photoe- 
mission experiments on doped materials show a d-wave 
anisotropy of the pseudogap in the superconducting state 
0j . In underdoped materials even in the normal state this 
pseudogap persists ||[|, wmcn is believed to cause the 
unusual non-Fermi-liquid behavior in the normal state. 
This emphasizes the importance of achieving a better un- 
derstanding of the superconducting phase, i.e. the phys- 
ical origin of the pairing mechanism, the nature of the 
pairing state and the character of low energy excitations. 

On a phenomenological basis the d-wave normal state 
pseudogap as well as the transition to a superconduct- 
ing state with a d-wave order parameter has been de- 
scribed within theories where short-ranged antiferromag- 
netic spin fluctuations mediate pairing in the cuprates 
11 • 

On a microscopic level it is believed that the Hubbard 
model or closely related models like the t-J model should 
capture the essential physics of the high-T c cuprates . 
However, despite years of intensive studies, these models 
remain unsolved except in one or infinite dimensions. 

Finite size QMC calculations for the doped 2D Hub- 
bard model in the intermediate coupling regime with 
Coulomb repulsion U less than or equal to the bandwidth 
W, support the idea of a spin fluctuation driven inter- 
action mediating d-wave superconductivity H. But the 
fermion sign problem limits these calculations to temper- 
atures too high to observe a possible Kostcrlitz-Thouless 



transition for the 2D system Q|. Another problem en- 
countered in QMC calculations is their finite size char- 
acter, which makes statements for the thermodynamic 
limit dependent on a scaling ansatz. 

These limitations do not apply to approximate many 
particle methods like the Fluctuation Exchange Approx- 
imation (FLEX) HI). Results of FLEX calculations for 
the Hubbard model are in agreement with QMC results, 
i.e. they show evidence for a superconducting state with 
d-wave order parameter at moderate doping for suffi- 
ciently low temperatures (|||. But the FLEX method 
as an approximation based on a perturbative expansion 
in U breaks down in the strong coupling regime U > W, 
where W is the bare bandwidth. On the other hand 
it is believed that a proper description of the high-T c 
cuprates in terms of the one-band Hubbard model re- 
quires U > W. 

Calculations within the Dynamical Mean Field Ap- 
proximation (DMFA) can be performed in the strong 
coupling regime and take place in the thermodynamic 
limit. But the lack of non-local correlations inhibits a 
transition to a state with a non-local (d-wave) order pa- 
rameter. The recently developed Dynamical Cluster Ap- 
proximation (DCA) |Ti]-[l| is a fully causal approach 
which systematically incorporates non-local corrections 
to the DMFA by mapping the lattice problem onto an 
embedded periodic cluster of size N c . For N c = 1 the 
DCA is equivalent to the DMFA and by increasing the 
cluster size N c the dynamic correlation length can be 
gradually increased while the DCA solution remains in 
the thermodynamic limit. 

Using a Nambu-Gorkov representation of the DCA we 
observe a transition to a superconducting phase in doped 
systems at sufficiently low temperatures. This occurs in 
the intermediate to strong coupling regime U > W and 
the corresponding order parameter has d-wave symmetry. 

Method A detailed discussion of the DCA formalism 
was given in previous publications |ll]-[l3| where it was 
shown to systematically restore momentum conservation 
at internal diagrammatic vertices which is relinquished 
by the DMFA. However, the DCA also has a simple 
physical interpretation based on the observation that the 
self energy is only weakly momentum dependent for sys- 
tems where the dynamical intersite correlations have only 
short spatial range. The corresponding self-energy is a 
functional of the interaction U and the Green function 
propagators. The latter may be calculated on a coarse 
grid of N c = L D selected K-points only, where L is the 
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linear dimension of the cluster of K-points. According 
to Nyquist's sampling theorem jL4j, this sampling of the 
reciprocal space at intervals of Ak — 2n/L implies that 
the DCA incorporates nonlocal dynamical correlations 
with a spatial range < L/2 and cuts off longer ranged 
dynamical correlations. Knowledge of the momentum 
dependence on a finer grid may be discarded to reduce 
the complexity of the problem. To this end the first Bril- 
louin zone is divided into N c cells of size (2ir/L) D around 
the cluster momenta K (see Fig|l]) . The Green functions 
used to form the self-energy E(K,w) are coarse grained, 
or averaged over the momenta K + k surrounding the 
cluster momentum points K (cf. Fig|I]). 
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FIG. 1. Choice of the N c = 4 cluster K-points (filled cir- 
cles), corresponding coarse graining cells (shown by different 
fill patterns) and a sketch of the d-wave symmetry of the order 
parameter. 

Thus, the coarse grained Green function is 

d(K,^H^£G(K + k,c) , (I) 

k 

where the sum runs over all vectors k = K + k within 
a cell around the cluster momentum K. Note that the 
choice of the coarse grained Green function has two well 
defined limits: For N c — 1 the sum over k runs over the 
entire Brillouin zone, G is the local Green function, thus 
the DMFA algorithm is recovered. For N c = oo the k- 
summation vanishes and the DCA becomes equivalent to 
the exact solution. The dressed lattice Green function 
takes the form 

G(k,L0) = (L0l-6 k T 3 ~£(K,Uj)y 1 , (2) 

with the self-energy E(k, to) approximated by the clus- 
ter self-energy E(K,w). To allow for a possible transi- 
tion to the superconducting state we utilized the Nambu- 
Gorkov matrix representation |T^] in ^ where the self- 
energy matrix E is most generally written as an expan- 
sion E = J2 . YiiTi in terms of the Pauli matrices t^. The 



diagonal components of E represent quasiparticle renor- 
malizations, whereas the offdiagonal parts are nonzero in 
the superconducting state only. 

Since the self-energy E(K,w) does not depend on the 
integration variable k, we can write 

G(K, u) = (wl - e K r 3 - E(K, w) - f (K, (3) 

where ck = N c /N e K+k- This has the form of the 
Green function of a cluster model with periodic bound- 
ary conditions coupled to a dynamic host described by 
T(K, u>). Here we employ the NCA to calculate the clus- 
ter Green function and self-energy respectively. A de- 
tailed discussion of the NCA-algorithm applied to the 
cluster model for the paramagnetic state was given in 
p2| . The NCA for the superconducting state has to be 
extended in order to account for the hybridization to the 
anomalous host, which couples cluster states with differ- 
ent particle numbers. 

The self-consistent iteration is initialized by calculating 
the coarse grained average G(K) (Eq. |f|) and with Eq. 
U the host function f (K), which is used as input for the 
NCA. The NCA result for the cluster self-energy E(K) 
is then used to calculate a new estimate for the coarse 
grained average G(K) (Eq. |l|). The procedure continues 
until the self-energy converges to the desired accuracy. 

Results We investigate the single particle properties 
of the doped 2D Hubbard Model 

H = H *« 4t c ja + Uj2 n nnn , (4) 
%j,a i 

where c\ (cj creates (destroys) an electron at site i with 
spin a and U is the on-site Coulomb repulsion. For the 
Fourier transform of the hopping integral tij we use 

= e — fi — 2i(cos k x + cos k y ) — At' cos k x cos k y , (5) 

accounting for both, nearest neighbor hopping t and next 
nearest neighbor hopping t'. We set t — 0.25eV and 
U = 3eV, well above the bandwidth W = 8t = 2eV. For 
this choice of parameters the system is a Mott-Hubbard 
insulator at half filling as required for a proper descrip- 
tion of the high-T c cuprates. 

To allow for symmetry breaking we start the iteration 
procedure with finite offdiagonal parts of the self-energy 
matrix E. As we mentioned above one expects the order 
parameter of a possible superconducting phase to have 
rf-wave symmetry. Therefore we work with a 2x2-cluster 
(N c = 4), the smallest cluster size incorporating nearest 
neighbor correlations. For the set of cluster points we 
choose K a ; = /tt, where I = 0, 1 and a = x or y. Fig.|l| 
illustrates this choice of K-points along with a sketch 
of the d-wave order parameter and the coarse graining 
cells. Obviously, for symmetry reasons, in the case of 
d-wave superconductivity, we expect the coarse grained 
anomalous Green function to vanish at the zone center 
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and the point (tt, tt). Whereas the anomalous parts at 
the points (0, tt) and (jr, 0) should be finite with opposite 
signs. 
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FIG. 2. (a) The local density of states (DOS) near the 
Fermi energy and the anomalous coarse grained Green func- 
tions at the cluster points (b) K = (0,0) and K = (n,n), 
(c) K = (tt, 0) and (d) K = (0, n) in the superconducting 
state. The nearest neighbor hopping integral t = 0.25eV, 
next nearest neighbor hopping integral t' = 0, bandwidth 
W — 2eV, the on-site Coulomb repulsion U = 3eV, tempera- 
ture T = 137K and the doping S — 0.19. The anomalous parts 
of the Green function (b)-(d) are consistent with a d-wave or- 
der parameter. 

Fig.|| shows a typical result for the local density of 
states (DOS) in the superconducting state along with the 
anomalous coarse grained Green function Gi2(K,w) = 

^/■WEk(( c K+it! c -(K+k)|}}w at tlie c l uster K-points 
for t' = 0, temperature T = 137K and doping 8 = 0.19. 
The anomalous coarse grained Green function vanishes 
at the cluster points (0, 0) and (tt, tt) but is finite at the 
points (tt, 0) and (0, it), consistent with a d-wave order 
parameter. Note that this result is independent of the 
initialization of the self-energy, i.e. an additional initial 
s-wave contribution vanishes in the course of the itera- 
tion. Thus a possible s-wave contribution to the order 
parameter can be ruled out. 

The finite pair amplitude is also reflected in the lo- 
cal density of states (DOS) depicted in Fig.||a, where we 
show the lower sub-band of the full spectrum near the 
Fermi energy. It displays a pseudogap at zero frequency 
as expected for a d-wave order parameter. 

Fig.| shows the DOS near the Fermi energy for the 
same parameters as in Fig.^], fixed temperature T — 
137K, but for various dopings. Obviously, the pseudo- 
gap size, measured as the peak to peak distance, as well 
as the density of states at the Fermi energy do not de- 
pend strongly upon doping. However the drop in the 



density of states from the gap edge to the uo = value 
first increases, reaches a maximum at about 19% doping, 
then decreases again. 
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FIG. 3. Density of states in a narrow region at the Fermi 
energy for the same parameters as in Fig.Q but for various 
dopings. The gap size and the density of states at uj = 
are independent of doping. Inset: Equal time coarse grained 
anomalous Green function Gi2(K,r = 0) at K = (ir,0). 

This behavior originates in the doping dependence of 
the anomalous Green function. In the inset we plot 
the coarse grained anomalous equal time Green func- 
tion Gi 2 (K,r = 0) = AV^Ek(c K +kT c -(K+k)i) for 
K = (tt, 0). This number as a measure of the super- 
conducting gap shows exactly the same behavior as the 
pseudogap in the density of states. 

The anomalous components Gi2(K, u>) and hence the 
pseudogap in the DOS become smaller with increasing 
temperature and eventually vanish at a critical temper- 
ature T c depending on the set of parameters. The phase 
diagram is shown in Fig.^. As a function of doping, T C (S) 
has a maximum T™ ax w 150K at 8 w 19% and strongly 
decreases with decreasing or increasing 6. The qualitative 
behavior of T c (6) in the calculated T—8 region agrees well 
with the generic phase diagram of the high-T c cuprates. 
Unfortunately, due to the break-down of the NC A at very 
low temperatures we are not able to extend the phase di- 
agram beyond the region shown in Fig^. This means 
in particular that we cannot predict reliable values for 
5 C (T = 0), beyond which superconductivity vanishes. 

The inset of Fig.|] shows the transition temperature 
dependence T c (t',S = const.) on the next nearest neigh- 
bor hopping amplitude t' for fixed doping 8 = 0.18. As 
compared to t' — T c strongly decreases with growing 
negative t' but increases for t' > 0. The shape of the 
phase diagram as well as the t'-dependence of T c can be 
qualitatively understood in terms of the phenomenologi- 
cal picture, where spin fluctuations mediate the electron- 
electron interaction, which then is strong at the antifer- 
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romagnetic wave vector Q = (jr, tt). 
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FIG. 4. Temperature-doping phase diagram for the 2D 
Hubbard model via DCA for a N c = 4 cluster. The near- 
est neighbor hopping t = 0.25eV, next nearest neighbor hop- 
ping t' = and the Coulomb repulsion U = 3eV. The error 
bars result from the finite resolution in temperature. Inset: 
Transition temperature Tc(t') for fixed doping S = 0.18 as a 
function of the next nearest neighbor hopping amplitude t' . 

In Fig.|5| we display the coarse grained spectra 
— — 3JmGii(K, uu) at K = (tt, 0) in the normal state 
(T = 290K) for the next nearest neighbor hopping t' = 
and t' = — 0.05eV (left and right hand side) for different 
dopings S. At the bottom we show the Fermi surfaces of 
the corresponding noninteracting systems (U = 0) in the 
first quadrant of the BZ. The diagonal thick solid line 
indicates the set of k-points which fulfill the nesting con- 
dition, i.e. which can be connected by Q to equivalent 
k-points in the opposite quadrant of the BZ. 

The doping dependence of the t' = and t' = —0.05 
spectra is qualitatively different. In the t' = —0.05 case 
the parts of the Fermi surface near K = (tt, 0) and (0, tt) 
fulfill the nesting condition roughly for the whole doping 
range, the quasiparticles couple strongly to the spin fluc- 
tuations and hence the corresponding spectra display a 
pseudogap at zero frequency over the entire doping range. 
The t' = spectra in contrast exhibit the pseudogap in 
the underdoped regime only (<5 — 0.05), where the spin 
fluctuations are strong, but show a quasiparticle peak at 
optimal doping 6 = 0.19, where the points on the Fermi 
surface near K = (tt, 0) and (0, tt) are far from being 
nested. The suppression of the density of states at the 
Fermi energy results in a suppression of superconduc- 
tivity and hence the transition temperatures drop with 
decreasing doping as well as decreasing t' < 0. For pos- 
itive t' we obtain similar spectra and Fermi surfaces as 
for t' = 0, but with a slightly enhanced density of states 
at the Fermi energy, resulting in higher transition tem- 
peratures. 



FIG. 5. Coarse grained spectral density — 9mGn(K,u) 
at K = (fT, 0) in the normal state (T = 290K) for the next 
nearest neighbor hopping t' = (left) and t' = — 0.05eV 
(right) for different dopings 5. The nearest neighbor hopping 
t — 0.25eV and the Coulomb repulsion U = 3eV. Bottom: 
Fermi surfaces of the corresponding noninteracting (U = 0) 
systems in the first quadrant of the BZ. The thick diagonal 
line indicates the set of k-points which fulfill the nesting con- 
dition. 

Summary We have used the recently developed DCA 
to study the long open question of whether the 2D Hub- 
bard model shows instabilities towards a superconducting 
state in the intermediate to strong coupling regime. We 
find conclusive evidence that at moderate doping a tran- 
sition to a state with offdiagonal long range order occurs 
and that the corresponding order parameter has pure d- 
wave symmetry. The corresponding temperature-doping 
phase diagram agrees qualitatively with the generic high- 
T c phase diagram. 
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